Prompt Quark Production by exploding Sphalerons 
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Following recent works on production and subsequent explosive decay of QCD sphaleron-like clus- 
ters, we discuss the mechanism of quark pair production in this process. We first show how the 
gauge field explosive solution of Luscher and Schechter can be achieved by non-central conformal 
mapping from the 0(4)-symmetric solution. Our main result is a new solution to the Dirac equation 
in real time in this configuration, obtained by the same inversion of the fermion 0(4) zero mode. 
It explicitly shows how the quark acceleration occurs, starting from the spherically 0(3) symmetric 
zero energy chiral quark state to the final spectrum of non-zero energies. The sphaleron-like clusters 
with any Chern-Simons number always produce Nf LR quarks, and the antisphaleron-like clusters 
the chirality opposite. The result are relevant for hadron-hadron and nucleus-nucleus collisions at 
large y/s, wherein such clusters can be produced. 



I. INTRODUCTION 

Although it would hardly ever be possible to observe 
instanton-induced events in electroweak theory, intense 
studies of such processes have started in this context . 
Even prior to these studies, the decay of the sphaleron so- 
lution |H has been studied in great details numerically j| . 
It was found that the Chern-Simons number at the exit 
state is frozen at some non-integer value. The question 
then is what happens to the relationship between the 
anomaly and the fermion chirality flow? Naively, the 
two concepts are related, yet the fermion spectral flow 
can only accomodate integer numbers. Should there be a 
generalization of index theorem for non-zero fields radi- 
ated to infinity? This question remained unanswered in 
general, although several points of view and partial re- 
sults have been expressed in the litterature. In this paper 
we would not address these questions in general, working 
out the solution explicitly, by finding exact solutions to 
the Dirac equation in exploding background fields, which 
demonstrates in detail how quark production and subse- 
quent acceleration occurs. 

Recently the whole discussion of such questions has 
shifted into the QCD domain. In refs it was 

suggested that nonperturbative configurations composed 
of instanton/antiinstanton play an important role in 
parton-parton scattering amplitudes at high-energy, and 
may account for some part of the soft pomeron: the esti- 
mates resulted in the slope and intercept in approximate 
agreement with data. Detailed discussion of a specific 
gluonic clusters which are produced have been made , 
with some proposed applications to high energy heavy ion 
collisions fig] . 

In the calculation of the cross section both the action 
of the gauge field corresponding to emitted gluons, or 
instanton-instanton interaction, can be naturally com- 
bined with the overlap factor with the incoming high en- 
ergy partons in the so called Landau method, in which 
the semiclassical approximation with singular gauge con- 
figurations is used . Pair singularities of the gauge 
field on the Euclidean time axis play the role of coordi- 
nate infinity in the original Landau work on quantum me- 



chanics. Those singularities are key to interpolating be- 
tween a vacuum configuration with zero energy, and the 
escape configuration with finite energy Q(T), where T 
measures the finite tunneling time. In a recent paper j^] 
we have shown that the gluon multiplicity associated to 
the escaping configurations at finite Q follows from the 
configuration at the sphaleron point Ms through a per- 
tinent rescaling of the sphaleron size and energy density. 

In the present work we continue along the same line, fo- 
cusing now on the quark production enforced by such ex- 
ploding sphalerons configurations. In section 2, we show 
how a conformal mapping from the 0(4) conformally 
symmetric solutions of the Euclidean Yang-Mills equa- 
tions leads to the 0(3) symmetric solutions by Luscher 
and Schechter jio) which describes the sphaleron explo- 
sion [ pT[ . In section 3, we construct the normalizable 
zero mode solutions to the 0(4) Yang-Mills background, 
which includes the instanton zero modes as a special case. 
In the central sections 4,5 we show that the solution 
to the massless Dirac equation in the Minkowski back- 
ground field of exploding sphalerons can also be obtained 
by the same conformal mapping, from the 0(4) Euclidean 
zero modes. We explicitly construct these states and 
show that at the initial Minkowski time t=0 those are 
zero energy states, while at asymptotically large time 
they reduce to a free quark or free antiquark of specific 
chirality. We also calculate the spectrum of the produced 
fcrmions. In section 6, we argue that pair production in 
the sphaleron (antisphaleron) background is through pair 
level crossing at t — 0. Our conclusions are in section 7. 



II. LS SOLUTION FROM 0(4) BY INVERSION 

In this section we show that the spherically symmet- 
ric 0(3) Luscher-Schechter (LS) solution describing 
the explosion of the sphaleron-type clusters JllJ can be 
obtained from the 0(4) symmetric solution through a 
pertinent inversion that breaks 0(4). Similar considera- 
tions for the 4 model were originally discussed in Jl2| . 
In our case, consider the 0(4) symmetric ansatz for the 
SU(2) gauge configuration 
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where 



Aapiv) = -Vanud„F(y) 



r Z(y) 

F(y) = 2 dffiO 
Jo 



ry\ The gauge configuration at t = satisfies Aq = and 
A1 = which corresponds to zero chromoelectric field 
(zero momentum) and finite chromomagnctic field. The 
t = point is a classical turning (or escape) point for 
the gauge configuration, where Euclidean and Minkowski 

(2) parts of the path join together, see more on that in [pd]| . 



with £(y) = \n(y/p). The 0(4) ansatz solves the Yang- 
Mills equations provided that 



III. THE 0(4) SYMMETRIC ZERO MODES 
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The solution with the turning point at £ — are express- 
ible in terms of elliptic sine and cosine, 
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for 2e < 1, while for 2e > 1 



f(0 = 5(1 
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In this section we will derive the generic zero mode so- 
lution to the Dirac equation in the 0(4) background (jl]). 
Our discussion in this section is limited to just one mass- 
less quark flavor, and will paralell the discussion in (ilf | 
for the instanton. For that, we introduce the conventions 



which arc related to the t'Hooft symbols througf 
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(ii) 



with Q/M s = 2e and 2m = 1 + 1/V2e. For Q = we 
recover the antiinstanton solution, while for Q = Ms we 
have the antisphaleron solution with 



m 2 ( cos(V20 



- 1 



(0) 



The instanton and sphaleron configurations follow from 
a similar construction with a dual background 77 — >• 77. 

The 0(4) solutions to the Yang-Mills equations with a 
turning point at £ = relates to the LS solution ]T3] with 
a turning point at t — by the conformal transformation 



(x + a) A 



2p 2 



(y + a)'' 



(y + <*)a 
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with a = (p, 0), which maps the sphere y 2 = p 2 onto the 
upper-half of the x-space. Indeed, under the conformal 
transformation (0), (Q) transforms as 



A a fi{x) 



(y + a) 2 
2p 2 



with y solving (0). The LS solution in Minkowski space 
follows by an analytical continuation X4 — > it, 

A a = f(0 Bptrn a 

g ((t - ip) 2 - r 2 )((t + ip) 2 -r 2 ) 
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4p(6 ai (t 2 -r 2 + p 2 ) + 2pe aij rij + 2n a rij) 
((« - ip) 2 - r 2 )((t + ip) 2 - r 2 ) 
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In terms of ( |lOf]ll| ) the Dirac operator in the background 
(0) reads 
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while in the dual background it reads 
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where the subsreipts c, s refer to the color, spin SU(2) 
matrices. The coupling g drops in the spectrum of D± 
making the quark states of order g°, while the gauge 
fields are of order 1/g. In the strict semiclassical ap- 
proximation, the quark effects on the semiclassical gauge 
configurations can be ignored. We note that the combi- 
nation 



P - 1 - - 1 + + 



(14) 



projects onto a color-spin singlet. In particular (U+ = 
(J2)aa is an eigenstate of P+, i.e. P+E/+ = U+. Finally, 
let ± L,R chiralities for which 75 ^f± = ±*&±. In the rj- 
field (D + ) or 77-ficld (D_) the zero modes carry specific 
chirality D±\I> ± =0. 

To construct the normalizable zero mode state in the 
configuration (|lj) , we note that only the positive chirality 
state is 'bound' 



(tr+ M + 2P+ a- c d, t F) 







(15) 
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while the negative chirality state remains free. The gen- 
eral solution to (U5f) can be sought in the form 



(16) 



up to an overall rigid color rotation and normalization 
(see below). Inserting ( |l6| ) into ( |15| ) and using the iden- 
tity 



it follows that 

□ $ = -2(^F)(a M $) 
which is solved by 

dy y 3 



(17) 
(18) 

(19) 



where F(£) is given by (|^). Using the identity (JTt]) for 
the norm of ( Jlq ) it follows that the overall constant is 
fixed by 



c= UJ A* ,2 (y) 



-1/2 



(20) 



The expression (16) thus defines the Euclidean normal- 
ized positive chirality (R) quark state in the 0(4) back- 
ground (|l|). Since the 0(4) Dirac spectrum is charge self- 
conjugate, the 0(4) background admits also a normalized 
negative chirality (L) anti-quark state in the 0(4) back- 
ground, i.e. V&L. The 0(4) background admits a pair of 
zero modes with opposite chirality. 



IV. SOLVING THE DIRAC EQUATION BY 
INVERSION 

In this section we use the conformal transformation (0) 
to map the pair of 0(4) zero mode states with turning 
point at £ = 0, onto a pair of 0(3)-symmetric solutions 
which solves the Dirac equation in the background LS 
configuration described in section 2. They are localized 
(normalizable) in 3-space at all times. They describe 
propagation of the produced quarks in Minkowski space- 
time into asymptotically free outgoing chiral quarks, with 
computable spectra. We will also show that those solu- 
tions in fact start with zero energy at the turning time 
t = 0: this implies that the sea contribution for positive 
and negative energy states is symmetric and cancels in 
the chiral charge. 

The solution itself is obtained by the following inver- 
sion formula 



Q+(x) = 74 fji rf * 
l/(y + a,y 



(y) 



(21) 



and it solves the (74 x ) Dirac equation in the (Euclidean) 
LS gauge configuration (||). Under the conformal trans- 
formation (Q) followed by 74 multiplication, the right- 
handed zero mode state (|16| ) remain right-handed, Un- 
winding y in terms of x and analytically continuing to 
Minkowski space yield 

Q+(*,r) = 



{{t + ipf - r 2 ) 2 
x (0 - it) + iS ■ x) e~ 2F ^ U+ , (22) 



with 



in 



(t - ip) 2 - r 2 
(t + ip) 2 - r 2 



(23) 



The normalization is fixed at t = by noting that F(0) 
0, i.e. 



Q+(0,r) = 



(p 2 + r 2 ) 2 
and demanding that in 3-space, 

1 
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(24) 



(25) 



which gives |C| = \/2p 3 /n. We note that the norm is 
independent of which LS configuration is chosen, making 
all mapped 0(3) states normalizable. 

The next question we address is whether the solution 
for the exploding sphaleron indeed starts with the zero 
energy state. Naively, at small t it does not look as a 
static solution at all. However, one can observe that the 
initial time derivative of (E2) fulfills 



-iQ + (0,r)= (3 - 8/ y 2 -^- f Q + (0 ) r ) . (26) 
This may be rewritten as a gauge phase 

-*Q+(0,r) = A(0,r)Q+(0,r) , (27) 
with SU(N) and U(l) gauge phases. The former is 



In A = —it- 



while the latter, Q + = e lt0 Q + , is 



e 



(3-8/(0))p 
p 2 + r 2 



(28) 



(29) 



Since Ao — at t — 0, the configurations (g2|) have 
zero energy modulo an SU(2)xU(l) gauge transforma- 
tion provided that /(0) is real. The SU(2) is related to 



Note that we are still using the notation C for the normal- 
ization of the zero energy state. 
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the color gauge group, and the U(l) to an additional lo- 
cal symmetry of the Dirac Hamiltonian H(t) = 747 • V 
at t = 0, namely 



M0 



H(t)e 



-ue 



H(t) - if 747 • ve 



(30) 



The solution thus looks static for the modified Hamil- 
tonian, and the difference between the transformed and 
original solution is absent at t = 0. We thus conclude 
that we do start from the zero energy state, masked by 
gauge transformations. All the 0(3) right-handed quarks 
at and above the sphaleron point are normalizable zero 
energy states to the Minkowski Dirac equation with the 
LS background. 




FIG. 1. The radial density distribution Q+(t,r) for 
t/p = 0, 1, 3, 6: the set of curves starts at the (anti)sphaleron 
point. 



At large times t ~ r ^> p with v = r — t, (|22j) simplifies 
-iC 



to 

Q+(*,r) 



At (v - ip)' A 



(1- B s 



2F(£)<?„-r 



with 



2 \v — ip 
We have used the relation 



-2F(S) 



(l-a s -f) = (l-a s -f) e 2F ^ 3 ^ 



(31) 



(32) 



(33) 



to display the exponential as an SU(2) hedgehog gauge 
rotation, since F{£) is purely imaginary asymptotically. 
The quark field is localized and normalized at t = and 
weaken asymptotically as 1 /£ up to a gauge transforma- 
tion. This is a hallmark of a radiation field which admits 



a normal mode decomposition as we will show below. 
The normalization is preserved by the time-dependent 
evolution since 



//■+00 
d 3 x|Q + | 2 (£,r)=7r|C| 2 y - 

The radial density is then given by 



+p 2 y 



= 1 



g+(t,r) =47rr 2 Tr(Q^Q + ) (f,r) 



that is 



/ x 16 

£+(<> r ) = — e 

7T 



4ReF(£) 



p 3 r 2 ((p 2 +t 2 ) +r 2 ) 
((t 2 +p 2 - r 2 ) 2 + 4r 2 p 2 y 



(34) 



(35) 



(36) 



at the (anti) sphaleron point as a function of the evolution 
time t/p — 0, 1, 3, 6. At the escape point t = 0, the zero 
energy state is fully localized in space. At large times 
t > p the energy density rapidly obtains a frozen form, 
and eventually the final quarks move luminally as free 
waves. A sample of density profiles is shown in Fig|l], 
where it is easily seen that the time needed to reach the 
asymptotic form is £* ~ p. 



V. SPECTRAL ANALYSIS OF THE OUTGOING 
QUARKS 

A normal mode decomposition allows for a spectral 
analysis of the final quark states released by the anti- 
sphaleron (sphaleron). Using the free field decomposition 
asymptotically, 



Q+(t,k) 



(27T)! 



(q«(£> 



-ikt 



ql(-fc)e 



+iki 



(37) 



allows the identification of the I +1 right-mover q# with 
the right-handed asymptotic chiral quark and the 1+1 
left-mover with the left-handed asymptotic chiral 
anti-quark. In the antispaleron background we have 



2k 



■Q+(t,k) 



(2tt) 3 / 2 

CVTfc 6(k) e - k P+ lkt (1 - a s ■ k) f/ 4 



so that 



q^(fc) = 2p 3 \l-e- k P{l-a s - k)U + 



hp , 



(38) 



(39) 



q/j follows from ( p9| ) by charge conjugation after mapping 
the 0(4) conjugate zero mode to the 0(3) zero energy 
state, 



q R (k) - 2p 3 



- k P{\ + a s -k)U^ 



(40) 



The (chiral) density of left antiquarks and right quarks 
are opposite asymptotically = — iir, with 
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Auk 2 
2k 



W R {k)\ 2 = p{2kpfe- 2k » 



(41) 



Fig. g shows the phase space distribution n R /p of the 
right quarks released at the sphaleron point as a func- 
tion of kp. The distribution integrates exactly to one 
produced quark, 



dk n^(fc) = 1 



(42) 



The quark spectrum is close to Planckian with an effec- 
tive temperature T — 2/p of about 300 MeV for a stan- 
dard p = 1/3 fm. Incidentally, this is close to the initial 
temperature of a quark-gluon plasma in the RHIC en- 
ergy domain. The released quarks carry asymptotically 
a total energy 



f°° 3 
Mj? = 2 / dkkn R (k) = - 

Jo P 



(43) 



which is small in the weak coupling limit under consid- 
eration, i.e. Mp/Mj = a/ V <C 1. We recall that in || 
the number of prompt gluons released was evaluated at 
1.1/ a s (p), which is parametrically large in this limit. For 
this reason, back reaction of fermions onto gluons can be 
neglected. 




FIG. 2. The spectrum of R quarks released at the 
Sphaleron point versus its momentum k in units of 1/p. 

The situation is however quite different for 'typical' 
QCD instantons and sphalerons, with p ~ 1/3 fm. In 
this case the produced 2Np = 6 quarks^] should be com- 
pared to ~ 3 gluons released classically from the explo- 
sion at the sphaleron point. The total energy carried 



3 We have ignored strange quark mass here, because the cor- 
responding scale at which the development occure is ~ 1/p ~ 
lGeV. 



by one flavor species would be 3/p ~ 1.8 GeV, while the 
whole available energy (the sphaleron mass) is only about 
3 GeV. Needless to say that such results follow from a 
wrong assumption: in this regime the fermions are not a 
small perturbation riding on the back of the gluonic wave. 
Their back reaction on the YM fields cannot be ignored. 
This imply a joint solution of the coupled Yang-Mills 
and Dirac equations is needed, better yet with vacuum 
polarization effects included. Although such analysis is 
still missing, we may speculate that the spectra would 
approximately look the same, with a rescalcd energy per 
quark appropriate for total energy conservation. 



VI. PAIR PRODUCTION BY SPECTRAL FLOW 

We have established that the release of an anti- 
sphaleron liberates N F LR pairs, while the release of a 
sphaleron liberates Np RL pairs. This mechanism oper- 
ates for all energies Q above and including the sphaleron 
point, and is signaled by a pair of energy levels crossing 
at t = (pair production) in the form of zero energy 
states, with their subsequent acceleration by pertinent 
gluoelectric fields. 

Compare this to what happens at low energies, when 
the gauge field configurations are well-separated instan- 
ton and anstiinstanton, and the production of 2Nf 
fermions is described by the well-known 't Hooft vertex. 
Although we have not shown it explicitly, we conjecture 
that the same number of fermions is produced at any 
energy, by continuity. 

On the other hand, the escaping gauge configurations 
above the sphaleron point were shown to carry fractional 
Chern-Simons number Q 



N(0) 



M s 



2/5 



(44) 



The non-integer character of ( f44| ) calls for a key question: 
How does the conventional Adler-Bell-Jackiw (ABJ) 
anomaly with integer fermion number work? The answer 
is as follows: The net chirality in the Dirac spectrum is 
carried by both the valence and sea parts, 



An(t) = (n fl -n L )(t) 

= Any(i) + An s (t) 



(45) 



At t=0 the Dirac spectrum associated to H(t) is C- 
conjugate in the 0(3) LS configuration except for the 
two zero energy states, 



An(0) = An 5 (0) = 2V F N(0) 



(46) 



For t > 0, one level dives in and one- level pulls out. The 
chirality count changes discretely by 2. In particular, 



Any(i) = 2N F 

An s {t) = 2N F {N(t) - 1) 



(47) 
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wherein the valence part is t-independent. Asymptoti- 
cally, the chiral polarization in the Dirac spectrum car- 
ried by the sea quarks vanishes Ans(oo) = 0, since 
the 0(3) asymptotic configuration is commensurate with 
plane waves, thus C-conjugate. The net result is 



An(oo) = Any(oo) 



2N f 



(48) 



which is the result explicitly obtained from the zero en- 
ergy states released. 

Some unsatisfied readers may ask how one can gener- 
alize index theorems, so that the number of level cross- 
ing would be calculable directly from the gauge field it- 
self, without the explicit solution of the Dirac eqn. Or, 
in a more practical form of the same, Does the confor- 
mal mapping exhausts all level crossings of the time- 
dependent Hamiltonian H(i) for all Q > Ms ? 

Interesting work toward answering these questions can 
be found in recent work where one can also find 
earlier references. In particular, it was found that the 
answer to the last one is negative. In a numerical anal- 
ysis presented there it was found that (in our units) for 
Q > 19.45 Ms other crossings in the LS background do 
occur. (However, these crossings are for practical pur- 
poses irrelevant since the cross section there is clearly 
too small.) 

To summarize: Sphaleron (antisphaleron) production 
leads to pair production of chiral quarks by spectral flow, 
and it produces one pair of quarks for each light flavor 
with unit probablity. This nonperturbative mechanism 
for prompt chiral production is not suppressed by powers 
of a s as in pQCD, or exponentially hampered by the con- 
stituent quark masses as in the string production mecha- 
nism. The qualitative effects of the current quark masses 
is to cause a power suppression in the production mech- 
anism as opposed to the exponential suppression in the 
string. 



VII. CONCLUSIONS 

We have started with zero mode solutions to the 0(4) 
configurations in Euclidean space that solves exactly the 
Yang-Mills equations. We then have shown that through 
a pertinent conformal mapping and analytical continua- 
tion they map exactly on the spherically 0(3)-symmetric 
LS solutions pi] with a turning point at t = describ- 
ing explosion of the YM sphaleron- like clusters jllj . We 
then have used the same conformal mapping to construct 
a class of normalizable solutions to the Dirac equation, 
in the LS background. We have shown that these solu- 
tions do correspond to the zero energy states at t = for 
all LS configurations with energies larger or equal to the 
sphaleron mass. 

By following their evolution in time, we have explic- 
itly shown that the zero energy states becomes at large 
time free Dirac outgoing waves with fixed chirality ±1 
and exactly 1 quantum. Thus sphaleron (antisphalerons) 



production liberates 2Np chiral quarks. We have iden- 
tified this liberation by a pair level crossing at t = in 
the Dirac spectrum. We have used the asymptotic zero 
energy states to construct the pertinent quark spectra 
produced. 

Many more things can be done based on the results 
reported here. One is an account of back reaction of 
fermions on the YM fields, as discussed in the previous 
section. Another is a projection of the outgoing waves 
to the wave function of outgoing hadrons, with many de- 
tailed predictions about exclusive and inclusive hadronic 
production to follow as we will report elsewhere. 

The present results are of interest to both hadron- 
hadron and nucleus-nucleus collisions at large y/s. Simple 
arguments indicate that while in hadron-hadron scatter- 
ing production via sphalerons, it is only a small part of 
the multiplicity, in nucleus-nucleus collisions it is not so 
and hundreds of sphalerons may be released, contributing 
substantially to the total entropy produced (T^| . In view 
of that, the production of 2Np chiral quarks per cluster 
may make the quark-gluon plasma produced in heavy-ion 
collisions quark rich, contrary to expectations based on 
perturbation theory or the color glass description Jig] . 
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